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We investigate the noisy Burgers equation (Kardar-Parisi-Zhang equation in 1+1 dimensions) 
using the dynamical renormalization group (to two-loop order) and mode-coupling techniques. The 
roughness and dynamic exponent are fixed by Galilean invariance and a fluctuation-dissipation the- 
orem. The fact that there are no singular two-loop contributions to the two-point vertex functions 
supports the mode-coupling approach, which can be understood as a self-consistent one-loop the- 
ory where vertex corrections are neglected. Therefore, the numerical solution of the mode coupling 
equations yields very accurate results for the scaling functions. In addition, finite-size effects can 
be studied. Furthermore, the results from exact Ward identities, as well as from second-order per- 
turbation theory permit the quantitative evaluation of the vertex corrections, and thus provide a 
quantitative test for the mode-coupling approach. It is found that the vertex corrections themselves 
are of the order one. Surprisingly, however, their effect on the correlation function is substantially 
smaller. 
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I. INTRODUCTION 

The Kardar-Parisi-Zhang (KPZ) equation represents 
one of the most prominent models describing nontriv- 
ial nonequilibrium dynamics This model equation 
constitutes one of the most thoroughly studied con- 
tinuum theories of kinetic roughening. It describes 
the height fluctuations /i(x, t) of a stochastically grown 
d-dimensional interface with a growth rate w(V/i) — 
A(V/i)^/2 depending nonlinearly on the local orientation 
of the surface, 



§ = -V^/^+^(V/^)^ + ^(x,<) 



(1.1) 



The {vV'^hyterm mimics a surface tension, and acts to 
smooth the interface, while the uncorrelated Langevin 
noise 77 (x, t) tends to roughen the interface and entails 
the stochastic nature of any growth process. Its first mo- 
ment vanishes, and its second moment is given by 

(77(x,t),7(x',0) =2i?5('^)(x-x')<5(i-t') ; (1.2) 

note that in general the coefficients v and D are not 
related in any simple manner, in contrast to near- 
equilibrium situations where Einstein relations connect 
damping constants and noise correlations. 



Dynamic scaling. The interface fluctuations are char- 
acteristically scale-invariant, i.e., the height profile ob- 
tained by a self-affine rescaling ft,'(x, i) — b~^h{bx,b^t) 
is, in a statistical sense, equivalent to /i(x, i). As a con- 
sequence, for sufficiently large Xq and to, such that the 
process is already beyond the initial transient region, the 
correlation function 

C(x, t) = ([/i(x + xo,t + to) - h{^o,to)f) (1.3) 

obeys the generalized homogeneity relation (x — |x|) 

C{x,t) = b~^^C{hx,hH) . (1.4) 

Upon choosing the scaling parameter 6 = 1/a; we obtain 
the dynamic scaling form 



C{x,t) = x^^C{t/x') 



(1.5) 



In the asymptotic limits t ~* and x 00, the scaling 
function C{t/x^) displays power law behavior and hence 



^ ' ' I Si^x/^ for a; ^ 



(1.6) 



The transverse wandering of the interface may be char- 
acterized by a perpendicular correlation length ^±{x) oc 
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•\/C(a;, <; = 0) cx with the roughness exponent x- The 
temporal increase of surface roughness is described by a 
parallel correlation length ^|| {t) cx t^/^ with the dynamic 
exponent z. 

Many growth phenomena show the above dynamic 
scaling of the interface fluctuations, but with values for 
the critical exponents different from those obtained for 
the KPZ equation. Nevertheless, the KPZ equation 
has become the starting point for our understanding of 
nonequilibrium dynamics and strong coupling behavior. 
Phenomenology of the KPZ equation. The phenomenol- 
ogy of the KPZ equation is now well known Be- 
low the lower critical dimension die — 2 there appear 
two renormalization-group (RG) fixed points, namely an 
infrared (IR)-unstable Gaussian fixed point and an IR- 
stable strong-coupling fixed point describing a smooth 
and a rough interface, respectively. For dimensions d > 2 
there exists a nonequilibrium phase transition from a 
weak-coupling phase for small effective coupling con- 
stants g — X^D/iy^, where the nonlinearity is irrelevant 
(in the RG sense), to a strong-coupling phase which 
seems to be inaccessible through perturbative meth- 
ods The scaling exponents in the strong-coupling 
phase have been determined by numerical methods 
and self-consistent mode-coupling approaches |||-||] . The 
results obtained from mode-coupling theory suggest the 
existence of an upper critical dimension due = 4 S . This 
result is supported by functional RG calculations Qjll) 
and renormalization group arguments [^;^. In the nu- 
merical simulations, however, the dynamic critical expo- 
nent z for the transient roughening of an initially flat 
interface is found to be smaller than zq — 2 for all di- 
mensions accessible to a numerical analysis i.e., there 
is no indication of any upper critical dimension. This dis- 
crepancy between mode-coupling theory and numerical 
results has yet to be resolved and constitutes one of the 
most important issues of current theoretical research. 
Mapping to other models. The KPZ equation is closely 
related to a variety of other problems ranging from fluid 
dynamics governed by the Burgers equation p3[ | to equi- 
librium systems with quenched disorder, namely directed 
polymers in random environments . Most of these 
mappings and relations are strictly valid for the one- 
dimensional case only. In order to assist the reader with 
the transfer of the results obtained in the main part of 
this paper to related systems, we provide a short account 
of some of the most important issues. 

The transformation v = — V/i leads to a Langevin 
equation for a randomly stirred fluid 



9v 
'dt 



A (v • V) V = i/V^v - V77(x, t) , 



(1.7) 



which in the case A = 1 represents a d-dimensional gen- 
eralization of the noisy Burgers equation p3[ . The long- 
time and large-distance behavior of the Burgers equa- 
tion, describing the dynamics of a vorticity-free velocity 
field, and the Navier-Stokes equation, characterizing an 



incompressible fluid, have been analyzed by Forster, Nel- 
son, and Stephen in the framework of dynamical renor- 
malization group theory to one-loop order [|l^. These 
authors have shown that the fluctuation-dissipation the- 
orem, valid in d = 1 only (see App. ^), together with a 
Ward identity resulting from the Galilean invariance of 
the fluid equation of motion allow the determination of 
the dynamic critical exponent z in d = 1 to be exactly 
z = 3/2. Their RG analysis has recently been extended 
to two-loop order ||l|,|,|l|l. 

Another model of surface roughening, which is gov- 
erned by the same nonlinearity as the KPZ equation, is 
the Kuramoto-Sivashinsky (KS) equation iQ. In con- 
trast to the KPZ equation the KS equation is completely 
deterministic 



dh 
'dt 



-vV'^h 



(1.8) 



and is characterized by a band of unstable modes at small 
wave vectors. (Note that v > 0.) Numerical simulations 
of the discretized one — dimensional KS equation have re- 
cently demonstrated that the large-scale dynamical cor- 
relations are described by the (1-1-1) -dimensional KPZ 
equation [Q. A derivation of the KPZ equation from 
the KS equation has also been given in Ref. jl^, where 
the effective parameters of the KPZ equation have been 
determined from the numerics of the microscopic chaotic 
dynamics of the KS equation. For d > 2, however, the 
results 20 1 are still controversial. 

Recently, Golubovic, and Wang succeeded in map- 
ping the equilibrium statistical mechanics of a two- 
dimensional smectic-A liquid crystal onto the nonequi- 
librium dynamics of the (l-l-l)-dimensional stochastic 
nonlinear KPZ (noisy Burgers) equation ||2^. Kashuba 
has shown that there exists a one-to-one relationship be- 
tween the Hamiltonian describing the nonlinear elasticity 
of a two-dimensional smectic-A liquid crystal and the 
Hamiltonian characterizing the long-range spin fluctua- 
tions in a two-dimensional planar ferromagnet subject 
to (two-dimensional) dipolar forces |^^. These relation- 
ships thus provide an interesting, exact approach to study 
the anomalous elasticity of smectic-A liquid crystals, as 
well as the spin fluctuations in the ordered phase of a 
dipolar planar ferromagnet in two dimensions, provided 
the corresponding KPZ growth model can be solved ex- 
actly, or at least to a high degree of accuracy. 

A number of somewhat more exotic relationships have 
been found very recently. E.g., the kinetics of the an- 
nihilation process A + B ^ with driven diffusion was 
mapped onto the (1-fl) -dimensional KPZ equation [ p5| , 
and the formal equivalence of the continuum limit of 
the Heisenberg equation of motion of a certain spin-1/2 
chain with the Fokker-Planck equation corresponding to 
the noisy Burgers equation was demonstrated [ p6[ . Be- 
sides these various mappings and relationships, which 
are valid in (1-1-1) dimensions only, the KPZ equation 
is also closely related to the dynamics of a sine-Gordon 
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chain the driven-diffusion equation p4| , and di- 

rected paths in a random media [|l4| . 
Invariances of the noisy Burgers equation. The one- 
dimensional KPZ equation is special in several ways. 
First, there is a huge list of mappings onto related mod- 
els as described above. Hence any advances in under- 
standing the growth model will have broad implications 
on many physical problems. Second, the noisy Burg- 
ers equation has two important "symmetry" properties, 
namely Galilean invariance and detailed balance. The 
Galilean invariance of the one-dimensional hydrody- 
namic equation (1.7), corresponds to an invariance of the 
stochastic growth model with respect to an infinitesimal 
tilt of the surface, /i^/i + v-x, x^x — Xvt. As a con- 
sequence of this symmetry, one finds that the amplitude 
of the nonlinearity A is invariant under RG transforma- 
tions, which in turn implies an exponent identity relating 
the roughness exponent x to the dynamic exponent z, 



X+z=2. 



(1.9) 



Whereas the latter invariance is valid for any dimension 
d, the detailed balance property of the KPZ equations 
holds in d = 1 only (see Appendix A) . It can be shown |^] 
that the Fokker-Planck equation corresponding to the 
(l-l-l)-dimensional KPZ equation has the stationary so- 
lution 



Vstih) oc exp 



'2D 



dx 



dh 
dx 



(1.10) 



this implies that the roughness exponent is x = 1/2, as if 
the nonlinearity we re e ntirely absent. Together with the 
exponent identity (O), one thus finds for the dynamic 
exponent z = 3/2. 

Scaling of the (l+l)-dimensional KPZ equation. As a 
consequence of the above inv arian ce properties of the 
nonlinear Langevin equation (1.1) one can show that 
the height-height correlation function obeys the follow- 
ing scaling law |M 



C{x,t) = Ax'^^F{\^/At/x'') 



(1.11) 



The argument of the scaling function is now dimension- 
less, and the scaling function itself is universal. It ac- 
quires the asymptotic form 



F(^\ _ / 1 for e - , 

^ I (e/2.g*)2^/^ for e -> oo 



(1.12) 



The RG fixed point of the (l+l)-dimensional KPZ equa- 
tion turns out to be a strong-coupling fixed point. As 
discussed above, despite this fact the roughness and 
the dynamic exponent are known exactly as a conse- 
quence of the particular invariance properties of the one- 
dimensional case. The scaling function F{^) has been 
calculated using a non-perturbative mode-coupling ap- 
proach . Striking agreement with the results of direct 
numerical simulations lES-ROl were found. 



The non-perturbative mode-coupling approach essen- 
tially consists in a resummation of the perturbation the- 
ory, such that all propagator renormalizations are prop- 
erly taken into account, while the vertex corrections are 
neglected. This is clearly a very ad-hoc and uncon- 
trolled procedure; nevertheless, mode-coupling theories 
have been remarkably successful in applications to many 
areas of condensed matter theory, such as structural glass 
transitions pl||, critical dynamics of magnets , bi- 

nary mixtures |3^j3^, and others |^^. In all those fields, 
it has been found that mode-coupling theory is capable 
of describing experiments in a quantitative manner. The 
factorisation approximation in the above mode coupling 
concepts is also known in the theory of hydrodynamic 
turbulence as Kraichnan's Direct Interaction Approxima- 
tion [||. 

The present work is motivated by this fact, and fur- 
thermore by the striking agreement of the mode-coupling 
results and those obtained from numerical simulations for 
the KPZ equation. In what follows, we will try to give 
a systematic analysis of the mode-coupling approach us- 
ing the field-theoretic formulation of Langevin dynam- 
ics [^^-Q. In particular, the fact that there are no 
singular two-loop contributions to the two-point vertex 
functions in perturbation theory in d = 1 strongly sup- 
ports the mode-coupling approach. As the IR singu- 
larities, i.e., the exponents z and Xi are exactly known, 
the self-consistent treatment is expected (and found) to 
reproduce the scaling functions to a high degree of accu- 
racy. In addition, we shall analyze vertex corrections in 
order to understand the range of validity of the mode- 
coupling approach. Our explicit results for the vertex 
corrections, as obtained from (exact) Ward identities, as 
well as from second-order perturbation theory allow for a 
quantitative estimate of the systematic errors enshrined 
in the mode-coupling approach. Since this specific type 
of self-consistent treatment is used in many areas of the- 
oretical physics, albeit under different nomenclature, we 
hope that this work will shed some light on its applica- 
bility, limitations, and possible extensions. 
Outline. The outline of the paper is as follows. In the 
subsequent section we summarize results from previous 
RG studies, discuss their relevance for the mode-coupling 
approach, and provide those explicit results which are 
needed in subsequent calculations. The formulation of 
the mode-coupling theory is discussed in Sec. Ill, as 
well as the solution of the self-consistent mode-coupling 
equations for the noisy Burgers equation. In addition to 
the scaling functions in the thermodynamic limit, finite- 
size corrections are explored. The size of the vertex 
corrections is estimated from the (exact) Ward identi- 
ties stemming from Galilean invariance, as well as from 
the explicit two-loop perturbational contributions. In 
the bulk of the present work, we shall refer solely to 
the (l-l-l)-dimensional KPZ equation; however, when- 
ever more general statements in d dimensions are pos- 
sible, this restriction to d = 1 is relaxed. We conclude 
with a brief summary and a discussion of some of the 
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remaining open problems. 



II. RESULTS FROM RENORMALIZATION 
GROUP THEORY 



r[h, h] - -F\j,j] + I d^x I dt{hj + h]) , 



where 



h = 5F/5j , and /i = 5F/5~j 



(2.4) 



(2.5) 



We start by reviewing some known results from pertur- 
bational renormalization group theory special- 
izing to 1+1 dimensions. This section also contains the 
explicit expressions for the vertex corrections to the two- 
point vertex functions to two-loop order. In this section, 
as well as in the Appendix, unrenormalized quantities are 
denoted by a subscript "0" . 



A. Dynamic functional 

We start with a brief description of the field-theoretical 
formulation of Langevin-type dynamics ]3^j37[ |. The 
stoc hast ic forces ?7(x, t) obeying (?7(x, i)) = and 
Eq. (1.2) can be taken to be Gaussian distributed, 



W[ri\ cx exp 



1 



4Z?o 



d'^x / dtrf{-K, t) 



(2.1) 



Using the equation of motion (1.1), we can eliminate 
the noise term; with an additional Gaussian transfor- 
mation introducing Martin-Siggia-Rose auxiliary fields 
h the ensuing probability distribution P[h] for the 
height fluctuations may be further linearized, and thus 
the original nonlinear stochastic equation of motion can 
be reformulated in terms of a generating functional 



dt 



jh + jh 



(2.2) 



with the Janssen-De Dominies functional given by 



J[h, h] 



d'^x dtlDo hh 



(2,3) 



Correlation and response functions can now be expressed 

as functional averages with weight exp Upon 

separating the dynamic functional into a quadratic and 
a nonlinear part, a standard perturbation theory can be 
formulated, where the cumulants Gj^ of the correlation 
and response functions are defined by functional deriva- 
tives of F[j,j] — In Z[j,j] with respect to the sources j 
and j, respectively. Vertex functions ^ are then ob- 
tained from the cumulants by a Legendre transformation. 



We finally note that the functional determinant originat- 
ing in the variable change from the noise fields 77 to the 
height fluctuations h serves to exactly cancel the acausal 
contributions to the perturbation series, thus leaving only 
those Feynman diagrams with correct time ordering in 
the response propagators ^^,0]. 



B. Two— point Vertex functions and renormalization 

We can now proceed to study the renormalization of 
the KPZ equation in one dimension. As discussed in 
detail in Ref. ||], the Ward identity stemming from the 
Galilean invariance of the Burgers equation shows that 
the nonlinearity A = Aq does not renormalize. This leaves 
the renormalization of the surface tension (diffusion coef- 
ficient) i^Q and of the noise correlation strength Dq, which 
may be inferred from studying the two-point vertex func- 
tions dq2Tj^^{q, w) and r^^(g, w), respectively; because of 
the fluctuation-dissipation theorem valid (only) in d = 1 
(see Appendix A), these coefficients are actually propor- 
tional to each other and must therefore renormalize in 
the same way. In Appendix B, we list the Feynman 
diagrams and the corresponding analytical expressions 
for Tj^f^{q,uj) to two-loop order (second-order pertur- 
bation theory in A), specializing the results of Ref. j|] 
to d = 1. Upon collecting these terms, splitting the 
vertex functions into regular and (UV) singular parts, 
^hh ~ ^^'^h ~'~ ^hh^ ^ eventually the following compara- 
tively simple results are obtained: 



+ ^ 1 



'piu! + VQq\ + i^oq'^ 



(2.6) 



(2.7) 



here we have introduced the abbreviations q± — (g/2)±p, 
q± = q±±k, and — dp/2TT. Note that the singular 
term stems entirely from the one-loop diagram (the ex- 
pression involving only one internal momentum p), while 
the (UV) singular two-loop contribution vanishes. The 
second-order term in the perturbation expansion thus 
yields merely regular corrections to the scaling functions. 
The second relevant vertex function can be written as 
F^^((7, w) — — 2£>oReF^^((7, w)/i'o'Z^, which allows us to 
define the wavenumber- and frequency-dependent diffu- 
sion coefficient as 
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iy{q,uj) 



2Dq 



confirming the validity of the fluctuation-dissipation the- 
orem of Appendix A [p[p^ . 

In evaluating those contributions which become singu- 
lar as the critical dimension die = 2 is approached, one 
has to be careful to choose a normalization point (NP) 
where either q or oj are finite, in order not to interfere 
with the IR singularities, which would also appear as 
poles in £ = d — 2 (for a more detailed discussion, see 
Refs. [|] and [Q). A convenient choice is NP: g = 0, 
iLol2v = K^; with go — X^Dq/vq one thus arrives at 



1 + 



90 



4 Jp k'^Z + p^ 



1+^ 



4 Jp k'^Z + p^ 



(2.9) 
(2.10) 



where Z is the renormalization factor for both v and D. 
The remaining singular integral is readily evaluated using 
the dimensional regularization scheme 



1 



d^l d/2 



(2.11) 



is a geometry factor. 



where Cd = r(2 - d/2)/2' 
and Ci — 1/2. Note that in this evaluation at fixed di- 
mension d — 1, no expansion with respect to e = d — 2 
was applied; the latter parameter was merely used to 
effectively count the singularities in the integrals that 
would appear at die = 2, when they are generalized to 
arbitrary dimension d. These ultraviolet (UV) poles may 
now be absorbed in renormalized quantities D = ZDq 
and ly = Zvq, with the renormalization constant 



Z = 1- 



8e 128e 
Defining the renormalized coupling 

50 

9 = 



(2.12) 



(2.13) 



Z^K ' 

we can now readily calculate Wilson's flow functions, 

C{g)=Kd^\oliiZ=-g/8, (2.14) 
(3ig)= Kd^log = g{d-2~ 2C) = .g(-l + .9/4) (2.15) 

in d = 1. Searching for zeros of the beta function yields 
the (IR) stable nontrivial fixed point 



5* =4, 

from which the critical exponents 

x= -Cig*) = 1/2 , 

z= 2 + C(<?*) = 3/2 



(2.16) 



(2.17) 
(2.18) 



can be deduced. Not e th at these explicit results fulfill the 
exponent sum rule (O); of course, as these exponents 
can already be determined from this identity and the 
additional constraint of the fluctuation-dissipation theo- 
rem (see Appendix A), this rather serves as a check for 
the calculations. Note that the remarkable cancellation 
of the singulartwo-loop contributions has been essential 
here from the diagrammatic point of view. 

In Ref . 1^ , the renormalization group approach is car- 
ried out in arbitrary space dimension < d < 4. For 
d > 2 an expansion with respect to e = d — 2 can be 
pursued, and was in fact recently carried through to ar- 
bitrary order in the perturbation series by Lassig . For 
d < 2, on the other hand, one may note that the fixed 
point coupling g* oc d approaches zero for d — > 0, and the 
results may be cast into an expansion about zero space 
dimension ifia. 



C. Two— loop scaling functions 

For later use, we now summarize the results from the 
second-order perturbation theory once more, albeit with 
some slight changes. First, we explicitly separate the 
zero- and one-loop contributions, and the two-loop con- 
tributions due to propagator and vertex renormaliza- 
tions. Second, we take "self-consistent" propagators, i.e., 
we generalize vq and Dq to a (/-dependent quantity ac- 
cording to Eq. (2_^), however neglecting its frequency 
dependence. This is in the spirit of the Lorentzian ap- 
proximation in mode-coupling theory, to be discussed 
below; its formal advantage is that the pole structure 
in the complex frequency plane remains unaltered, and 
therefore the results from Appendix B may be readily 
generalized. The zero- and one-loop contributions to 
r^^(g, 0) thus read (see Fig. ^,b): 



iih 



(9,0) 



1 



iq+)ql + u{q-)q^ 



(2.19) 



similarly, the two-loop contribution due to propagator 
renormalization (Fig. 0c-f) becomes 



q^_ 



v{q+)q\ + u{q+)q\ + v{q^)q^^ 



(2.20) 



while the result for the two-loop contribution due to ver- 
tex corrections (Fig. ^-j) is 

rg:^(<Z,0) ^ -q'\' I I , . , , , (2.21) 

Jph v{q+)qX + v{q-)q- 

{v{q+)q\ + v{q_)qi\{v{q+)q\ + v{q-)qi + i'(fc)fc2] ' 



5 



III. MODE COUPLING THEORY AND VERTEX 
CORRECTIONS 

In this section we study the mode coupling approxima- 
tion for the Burgers equation. For readers not familiar 
with the dynamical functional approach discussed in the 
previous section, we start by a derivation of the mode 
coupling equation using a perturbation theory for the 
equation of motion. 



A. Perturbation series and mode coupling equations 



In Fourier space with the equation of motion (|l.l| ) 
reads 



h(k, uj) = Go(k, w)ry(k, uj) + Go(k, oj)j{k, oj) 

+ iGo(k,c.) / V^^]^ hik+,u;+)hik.,u;_) 
^ Jq.u 



(3.1) 



where Go(k, ~ \/{v\i^ 
tor", K 



(0) 



ijj) is the "bare propaga- 
-Aki • k2 is the "bare vertex", and 
k± = k/2 ± q, u)± = ijj/2 ± /i. The noise r] is assumed to 
be Gaussian and uncorrelated, given by the weight func- 
tion Eq. (2.1). A small external perturbation j(k, lS) has 
also been included in Eq. (|3.lD, and will be used to gen- 
erate the response functions. Typically, the quantities 
of interest are the noise-averaged two-point correlation 
function 

(/i(K)/i(K')) - G(K)5(K + K') , (3.2) 
and the noise-averaged linear response function 



/ 5h{K) 



G(K)5(K-K') , 



(3.3) 



where (k, w) is abbreviated by K and (5(K + K') = 
{2TiY+'^5'^{k + k')5{u} + to'). These are special cases of 
the general Green's function 



Gi, 



i(— Pi; ■ • • ; ~Pr 
'(5'"(MKi)- 



|Ki; . . . ; K„) 
■h{Kn))' 



(3.4) 



where the subscript "c" denotes the connected part. 
In this notation, the two-point correlation function is 
{h{K)h{K')) = Go,2(-|K;K') and the linear response 

fimction is ((5/i(K)/(5j(K')) = Gi,i(-K'|K). From 
Eq. (B.4), it is clear that Gm,Q = 0. The above definition 
of the Green's functions is identical to the one used in 



the dynamic functional formalism (Sec. II A). 

One approach to study the Green's functions Gm,n 
perturbation theory. For F^"^ = 0, Eq. (3J) is just the 
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linea r di ffusion equation. For ^ 0, the solution of 
Eq. (3.1) may be obtained iteratively by a perturbation 
expansion in powers of V^^\ For example, the lowest 
order correction to the response function is 



Gi(K) = Go(K) + Go(K)Ei(K)Go(K) , (3.5) 
where Go(K) = 2D\GoiK)\'^ is the "bare correlator" , and 



Si(K) 



(0) 



Q 



k+;k 



Go(K_)Go(K+)V; 



(0) 
k+;k 



(3.6) 



is the one-loop renormalization of the "self-energy" . 
Similarly, the lowest order correction to the correlation 
function is 



Gi(K) = 2i?i(K)|Go(K)|^ 



(3.7) 



(0) 

k+;k_ 



(3.8) 



where 

D^{K) = D + \j^ vi%_C^{K^)C^{K+)V, 

is the one-loop renormalization of the "noise spectrum" . 
Unfortunately, such perturbation series diverge in the hy- 
drodynamic limit k, — > 0. One way to proceed is to 
perform a renormalization group analysis. It turns out, 
however, that there is no fixed point that can be obtained 
in a controlled e expansion with e = 2 — d below d — 2 
dimensions . Hence, a non-perturbative method is re- 
quired to treat the KPZ problem. One approximation 
which has been frequently used is to replace the bare 
propagator Gq and bare correlator Gq in Eqs. (3^) and 
( |3.8[ ) by the renormalized functions G and G while keep- 
ing the vertex V'^^^ unchanged. This is known as the 
mode-coupling approximation (or Kraichnan's Direct In- 
teraction Approximation), and it leads to the following 
closed set of integral equations, 



(0) 

k+;k ' 



I](K)=- / G(K_)G(K+)F, 
Jq 

DiK)=D + ^I^V^l]^_C{K.)CiK+) 



(0) 



k+;k_ ' 



(3.9) 



(3.10) 



where S and D are defined by G and G through 

G-i(K)=G„-i(K)-S(K), (3.11) 
G(K) = 2Z?(K)|G(K)|2 . (3.12) 

Of course, as this procedure neglects any vertex renor- 
malizations, it constitutes a partial sum of the perturba- 
tion series only, and as a-priori no information is avail- 
able about the size of the missing contributions, it clearly 
constitutes an uncontrolled approximation. Neverthe- 
less, the mode coupling theory has been quite successfully 
applied in many areas of condensed matter theory, as 
mentioned in the introduction. It was first applied to the 
KPZ problem by van Beijeren, Kutner, and Spohn to 
get the scaling exponents x and z in d — I. Recently, the 
mode-coupling equations were solved numerically to ob- 
tain the entire function G(k,a;) in 1 -t- 1 dimensions [ p7| , 
and striking agreement to the scaling function obtained 
by direct nu merica l simulations ||2^ were found (for de- 
tails see Sec. [II C below). This result is very surprising 
and prompted us to study the mode-coupling theory in 
more detail. In what follows, we will try to give a sys- 
tematic analysis of the mode-coupling approach using 
the field-theoretic formulation of Langevin dynamics. 
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B. Dynamic field theory, Ward identities, and vertex 
corrections 

The starting point of our study is the response func- 
tion, which can be formally obtained by differentiating 
Eq. (3.1) with respect to the perturbation j(K'). We 
obtain 



Gi,i(-K'|K) = Go(K)5(K-K') 



(3.13) 



Mode coupling theory amounts to expressing Gi^2 hi 
terms of the lower order functions Gi,i and Go, 2. 

In order to analyze the Green's functions systemati- 
cally, we turn to a functional integral method described 



in section ||. From the gen erat ing functional, Eq. ( |2.2| ) 
the Green's functions, Eq. (|3.4|) , can be easily obtained 
as the functional derivatives of = log Z[j,j]. 

By taking derivatives of Eqs. ( |2.5| ) and using Eq. (|]J), 
it is straightforward to relate the Green's functions Gm.n 
to the vertex function Tm,nS, e.g., 

ri,i(P|K) = ri,i(-P) = G-\P)S{K + P) . (3.14) 

The two-point correlation function can also be easily 
found. It has the form of Eq. (3.12), with 



r2,o(Pi, P2|-) = r2,o(Pi) = -2D{p,)S{Pi + P2) 



(3.15) 



All higher order Green's functions can be written as products of G(K), G(K), and the higher order vertex functions. 
For example. 



G2,i(Pi; P2IK) = -G(K) ri,2(K|Pi; P2) G(-Pi)G(-P2) , 
Gi,2(P|Ki;K2) = -G(Ki)G(K2)r2,i(Ki;K2|P)G(-P) 
-G(Ki)G(K2) ri,2(Ki|K2; P) G(-P) 
-G(Ki)G(K2)ri,2(K2|Ki;P)G(-P) . 



(3.16) 
(3.17) 



Using Eq. ( |3.17D in Eq. ( |3.13| ), we obtain 

Gi4(-K'|K)= G(K)(5(K - K') 

= Go(K)J(K-K')-Go(K) J V^l]^_G{K^) 



'Q 

G(K+)ri,2(K_|K+; -K') + iG(K+)r2,i(K_; K+| - K') 



G(K') 



(3.18) 



for the full response function. Note that it has the form 
G(K) = Go(K) + Go(K)S](K)G(K), where S(K) is the 
self energy defined in Eq. (3.11). If we write the vertex 
functions as 

ri,2(K|Ki;K2)=r„(Ki;K2)5(K-f K1+K2) , (3.19) 
r2,i(Pi; P2IK) = rb(P2; K)J(K + Pi -t- P2) , (3.20) 

then the self-energy becomes 

E(K) = - J^V^l\^_G{K^)C{K+)ViK+;K) , (3.21) 
where 

y(K+;K) = r„(K+;-K) 

G-i(-K- 



The additional term does not change S(K) in Eq. (3.21) 



r5(K+;-K) 



(3.22) 



4i?(K+) 

denotes the "renormalized vertex function" . It will be 
useful to write the vertex function in a slightly different 
form, 



F(K+;K)=rjK+;-K) 

G-HK+) + G-i(-K+) 



rb(K+;-K). (3.23) 



because its poles, from G(K_) and G(— K+), are on 
the same side of the complex frequency plane. Hence 
the frequency integral for t his addition al t erm yields 
zero. Comparing Eq. ( 3.21 ) with Eq. ( |3.9|) , we real- 
ize that the mode-coupling equation becomes exact if 



F(K+;K) ~ k- section III D| , we will show that 



this equality in fact does not hold. Yet, by exploiting 
a number of identities relating the different vertex func- 
tions, we shall show that the correction to y(K_|_;K) is 
small in the limit K — > 0. This is hopefully the first step 
in understanding the puzzle of why the mode coupling 
theory works so well, at least in the case of the noisy 
Burgers equation. 



C. Numerical solution of the mode coupling 
equations 



In this section we present the numerical solution of the 
(1-f l)-dimensional KPZ equation. In view of the results 
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from the preceding section it is convenient to define a 
generalized kinetic coefficient D(k,uj) and a generalized 
"surface tension" v{k^ijj) by 



1 



—iuj + v{k, to) ' 



C{k,Lu) 



2D{k,uj) 
uj'^ + [i^{k,uj)Y 



(3.24) 



(3.25) 



Then, the self-consistent equations for the correlation 
function C(k,uj) and the response function G{k,uj) in 
Fourier space are given by 



v{k,uj)^y klkk^C{k+,uj+)G(k^,uj^) , (3.26) 
4 , 



D{k,uj) 



klklC{k+,uj+)C{k^,uj_) . (3.27) 



For the numerical solution of the mode-coupling equa- 
tions it is much more convenient to study the intermedi- 
ate correlation and response functions, defined by full or 
half-sided Fourier transforms, respectively. 



(3.28) 
(3.29) 



C(fc,w) = / dte*"*C(fc,t) , 

•J —OQ 

G{k,uo) = / dte"^^e{t)G{k,t) 



Inserting into Eqs. ( 3.26 )-( 3.27 ) we get for the general- 
ized kinetic coefficients 



v{k,t) = AM kXkk^C{k+,t)G{k^,t) , 

Jq 

2 



D{k,t) 
where 



klk^_G{k+,t)C{k-,t) 



(3.30) 
(3.31) 

(3.32) 
(3.33) 



It is again important to realize that for the (1+1)- 
dimcnsional KPZ equation there exists a fluctuation- 
dissipation theorem (FDT) which relates the generalized 
kinetic coefficient D{k,t) with the generalized surfa ce 
tension h'{k, t). As shown in Appendix A (see Eq. (All)), 
the following identity holds 



/oo 
dte''^*D{k,t) , 
-OO 
nOO 

v{k,oj) = / dte"^^e{t)i^{k,t) 



G{k,t) ^ —e{t)G{k,t) 

This allows one to rewrite v^k, t) as 



v{k,t) 



2v 



G{k+,t)G{k^,t). 



(3.34) 



(3.35) 



Together with Eq. (3.24), which can be written as 



|-G(fc,t)==-/ dTiy{k,T)G{k,t-T), (3.36) 
c** Jo 

one now has a set of self-consistent equations for the ef- 
fective surface tension and the response function. 
Scaling analysis of the mode coupling equations. We look 
for the solutions of the scaling form 



iy{k,Lu) — Xvk^v{il)) , 
D{k,uj) = XDk-^'viiJj) , 



(3.37) 
(3.38) 



where we have defined the scaling variable Co — 
Lu I Xvk^ . The corresponding scaling forms for the 
Fourier-transformed quantities read 



iy{k,t) = {\v¥-fv{i) , 
D{k,t) = \^Dvk'--^'v{t) 



(3.39) 
(3.40) 



with the scaling variable t = Xvk^'t. For the response 
function the scaling analysis leads to 



G{k,uj) 
G(k,t) 



1 



Xi^k^ 
G{i) 



Giu) , 



(3.41) 
(3.42) 



Inserting the scaling forms, Eqs. (.3.37)-(3.40), into the 



mode-coupling equations implies for the dynamic expo- 
nent z — 3/2, and leads to the following self-consistency 
equations for the generalized kinetic coefficient and the 
response function: 



Hi) 



1 

2^ 



dxG{x+i)G{xJ) , 



^ G{i) = - f dTv{T)G{i-T) 



dt 



(3.43) 



(3.44) 



where x± = 1/2 ± a;, and the effective coupling constant 
is given by 



(3.45) 



Note that the amplitude A is arbitrary. We have chosen 
it to be equal to the effective coupling constant in order 
to simplify the scaled mode-coupling equations. 



It can be shown analytically from Eq. ( 3.43 ) that 
the scaling function for the generalized surface tension 
v{i) shows a power-law behavior v{t) = Di~'^^'^ with 
D Ri 0.1608 for small times t < 10^^. Since the response 
function is almost constant for small t, one finds from 



Eq. (3.44) that 



G(t) = exp • 



-Cgaussi^/'} for £< 10- 



(3.46) 
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with Cgauss = 9j7/4 « 0.3619. In Fig. | the numeri- 
cal solutions for the scaling functions z>(t) and G{t) are 
depicted, as well as the results from the Gaussian ap- 
proximation (3.46). 
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FIG. 1. Scaling functions for the generalized surface ten- 
sion vit) and response function G{t) vs. the scaling variable 
i — \{D /vY^^k^^'^t. The poi nt -da shed line represents the 
Gaussian approximation, Eq. (3.46), for the response func- 
tion, which is obtained from the analysis of the mode coupling 
equations at small times. 



Truncated correlation function in real space. Another 
quantity, which is easily accessible by numerical simula- 
tions, is the truncated correlation function in real space. 



C{x.t)= j |^2[l-e"-]C(t,t) 



(3.47) 



where F(0) = 1 since G(0) = 1. Conforming to the def- 
inition in Eq. (1.11) one gets A = D/v. The universal 
scaling function f(£) and is shown in Fig. 0. 
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FIG. 2. Scaling function F{^) for the truncated corre- 
lation function in real space versus the scali ng variable 
5 = XA^I'^t/x'^''^. The empirical forms, Eq. ( |3.50| ), are shown 
as the dashed curves. The dimensionless coupling constant 
can be read off from the crossover point of F{£^). 

The dimension less argument of F has the form de- 
manded by Eq. ( 1 . 1 1| ) with z — 3/2. The dimension- 
less coupling constant can be read off from the crossover 
point of F{£) (see Fig. |); we obtain g* = 0.87. This 
result can be checked more precisely in simulations by di- 
rectly looking at the scaling amplitudes. Our work thus 
predicts that if 



then 



C{x = 0,i) = 0.70(AA2t)2/3 



(3.48) 



(3.49) 



The numerical error is less than ±1%. The amplitude 
0.70 ± 0.1, extracted from the mode couphng equations, 
agrees rather well with results from numerical simula- 
tions, which find an amplitude of 0.712 ± 0.003 |^ and 
0.725 ± 0.005 [|9|, respectively. In Ref. ||| a empirical 
form for F{^) has been given, which fits the data from 
the numerical simulation quite well. We find that the 
mode-coupling result is also quite well approximated by 
the same empirical forms (dashed curves in Fig. |^) 



1 + 4.22 exp {-3. 82^3/2| 
0.7^2/3^0.43^'/' 



for ^> Co, 



(3.50) 



where ~ 2.5 (see the dashed curves in Fig. ||), but with 
somewhat different numerical values for the coefficients. 
Note that the dashed curves are almost indistinguishable 
from the solid line; in order to make the dashed curves 



visible we have plotted the asymptotic forms in Eq. 3.5C 
for values smaller and larger than respectively. In 



9 



summary, the mode coupling result for the scaling func- 
tion F(( ^ ) a grees with the results from numerical simula- 
tion pq , |2t^ within a few percent. 

Finite size effects. Note that the above results are valid 
for very large systems in the steady state. Transient be- 
haviors such as the growth of interfacial width starting 
from flat initial conditions may well be more compli- 
cated . They may also be computed using the mode- 
coupling theory (with Fourier-Laplace transform to in- 
corporate the initial conditions; however, the procedure 
becomes more cumbersome. 

Nevertheless, we can say something about the behav- 
ior of systems of finite size L already on the basis of 
our results for the steady state. In principle, the cor- 
relation function and response function are now explic- 
itly L-dependent. They may be described in terms of 
Di^{k,uj) ~ i'i^{k,uj) ^ L^^^f{Lj,kL), where / is now 
the solution of a two variable integral equation with the 
initial condition D^^aik^Lu) = D (the bare value) and 
similarly for i>. In this way, one would obtain the ex- 
plicit functional renormalization of various quantities as 
we look at larger length scales L. The flow behavior of 
D and v described by usual recursion relations are re- 
covered from the L dependence of Di(A; = 0, w = 0) and 
I'Lik = 0,uj = 0). The asymptotic form ^ vj^ ^ L^l'^ 
is of course the expected one given the exponents x ^ 
pO| . Here we want to emphasize that the self-consistent 
equations provide a connection between the microscopic 
and macroscopic (renormalized) theory. 

If the flow of these functions is already well advanced, 
i.e., for times much larger than the initial time t = 0, 
where the interface was absolutely flat, our results for the 
steady state can also be used to get approximate results 
for the "transient behavior" of a finite size system. Note 
that with "transient behavior" we are not referring to the 
transients starting out from an absolutely flat interface, 
but to transient behavior after some initial relaxation. 

The interface width in a system of finite size L is de- 
fined by 



wl{t)^l^h(x,t)-h{xM 



(3.51) 



Upon assuming that the spectrum of the height function 
is only slightly modified by finite size effects (and/or af- 
ter some initial transient), the interface width can be 
approximated in terms of the correlation function in the 
steady state 













oo 27r ^ 


4D 


/■+°° dk 1 


V 





27r/L 



dk 

2^ 



C{k,Lo) 



[1-G(fc,t)] 



(3.52) 



for periodic boundary conditions. Inserting the scaling 
laws Eqs. (3.37)-(3.40) for the correlation functions one 
finds 



wlit) = ^Lf{f) 



where f = AA/—i ^/^t and 



(3.53) 



f{r) 



dxx 



I" J27r 

Asymptotically one gets 



1 - G(f a;3/2) 



2t^\ _ } ClL 



for t ^ oo 
C^t^/^ fort^O . 



(3.54) 



(3.55) 



Just as in the discussion of the steady state correlation 
function one can define an universal amplitude ratio by 



R ^ 



Ct 



(3.56) 



and rewrite the scaling law in terms of this ratio 

WL{t)^CL^/lW{T) , (3.57) 

where r = XCLtL-^'^, Cl = f{oo)D/v with f{oo) « 
0.101. The scaling function W{t), shown in Fig. ^, has 
the following asymptotic behavior 



Miri^\ / 1 for r ^ oc 



(3.58) 



The ratio R is found to be i? = 3.8, which is in reasonable 
agreement with i?cxp ~ 3.45 ± 0.05 found in muncrical 
simulations BOtl. 



W(T) 



periodic boundary conditions 



10 




10" 



FIG. 3. Scaling function W{t) for the interface width in 
a system of finite size L for periodic boundary conditions 
and free boundary conditions, respectively. The dashed lines 
are approximation for small scaled times r = XChtL'^^^ , 
W{t) ~ 3.8r^'''^ for periodic boundary conditions and 
W{t) ~ 2.4r^/'^ for free boundary conditions. 
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If one uses free instead of periodic boundary condi- 
tions, o ne ha s to r eplac e the lower bound of the integral 
in Eqs. (B.52) and (3.54) by tt. The resulting scaling func- 
tion for the interface width is also shown in Fig. and 
we find, as already noted in Ref. that the asymp- 
totic behavior at small times is given by W{t) — R^t'^^'^ 
with Rf = i?/22/3 « 2.4. 



D. Vertex corrections 



limit h,h ^ 0, and recalling the definition of the vertex 
functions, we obtain the following Ward identity 



d 



iuj lim Trrri,2(qi,/ii I ki,wi;k,a;) 
= A<^ qiri,i(qi,^i +UJ I ki,wi) 



-kiri,i(qi,/zi I ki,uji +uj] 



(3.64) 



As we have seen in the preceding section mode cou- 
pling theory is equivalent to a self-consistent formulation 
of the perturbation series, where all propagator rcnor- 
malizations are taken into account, but vertex correc- 
tions have been neglected. Nevertheless, there is quite 
an excellent agreement of the mode coupling results with 
numerical simulations p8[ |. It seems that there is some 
hidden small parameter, which remains to be identified. 
In this section we address this problem and analyze the 
magnitude of the vertex corrections. 

It is known that the KPZ equation is invariant under 
a Galilean transformation of the form 



ft,'(x, t)= /i(x + Avf,t) + 
^'(x,i)= /i(x + Avi,t) , 



(3.59) 
(3.60) 



corresponding to an infinitesimal tilt v of the surface. 
This invariance leads to Ward identities, connecting the 
two- and three-point vertex functions ^J, which imply 
that the nonlinearity A is not renormalized, and that 
there is an exponent identity x + z = 2. 

Since the Ward identities relate the three-point with 
the two-point vertex functions one may hope that they 
also give some information on the magnitude of the ver- 
tex functions. Recently, it has been shown by Lebedev 
and L'vov |4^ that the KPZ equation is invariant under 
the generalized Galilean transformation 



/i'(x',t)=/i(x,t) + ^.x', 

h' {-k! ,t)^h{-K, t), 



(3.61) 
(3.62) 



with x' = x — AC and where C, is an arbitrary function of 
time but not of coordinates x. Since the generating func- 
tional for the vertex functions T[h, h] is invariant with re- 
spect to the above transformation, one finds the following 
Ward identity 



dt 



Ak-C 



5T 



Sh(k,t) 
ST dC d 



6T 



Sh(k,t) 



h(k,t) 



Sh{k, t) dt dk 



(5('*)(k) 



. 



(3.63) 



Taking functional derivatives of the above equation with 
respect to ft.(— qi, — /ii) and /i(— ki, — tJi), then taking the 



Similarly, by taking derivatives of Eq. (3.63) with respect 
to h{-qi, -^i) and /i(-q2, -^^2), we get 



d 



iuj lim — r2,i(qi,Ati;q2,At2 | k,u;) 

k^O OK 



= A<^ qir2.o(qi,/^i +'^;q2,Ai2 1 •) 



+q2r2,o(qi,Mi;q2,M2 + w | •) 



The general Ward identity reads 

icjVkr™,„+i({Qi}|{Ki}; k, w)|k=o 

m 

= X'^qjr^,n{{Qi} + ujej\{Ki}) 

n 

-FA^kjr„,„({Q, \{K,}+ujej) , 



(3.65) 



(3.66) 



where we have defined {Qi} = qi, /ii; . . . ; qm, Mm, 
{Ki} = ki,wi; . . . ;k„,w„, and {Q^} + we^ = 
qi, /ii; . . .]qj,fij + uj- . . . ; q^, Mm- I nsert ing the above 
Ward identities, Eq. (3.64) and Eq. (3.65), into the ex- 
pression for the vertex correction, we find 



VkV^(K+;K) |k=o= — q< -icj + j/*(q, /!_)-;/* (q, m+) 



'^(q, M+)+t^*(q,M+) 

2D{q,fi+) 



:i?(q,M-)"^(q,M+)] , (3.67) 



where fj,± = fj,± ^. The first term in the latter equation 
corresponds to the bare vertex, which is real. The correc- 
tions to this bare vertex result from the imaginary part 
of the frequency dependent surface tension v(k, lu) in the 
second and third term of Eq. (3.67). In addition, the 
renormalized vertex contains an imaginary part result- 
ing from the real part of the generalized surface tension 
and the noise amplitude D{k,Lu). 
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where ft = fi/Xiyq^^"^ and uj = uj/Xvq^/"^. As can be in- 
ferred from Fig. ^ the vertex corrections may be as large 
as 40% at certain values of the external frequencies. If 
we take the integral over all frequencies as a measure of 
the vertex correction, however, we find that it is only of 
the order of a few percent or even less. 



E. Vertex corrections from the two— loop 
contributions 



In this subsection we study the vertex corrections 
resulting from two-loop diagrams in Lorentzian ap- 
proximation. With the ansatz ^{q) — I'iqTLo = 
0) = t'ior(A/\/27r)g^^^ the mode-coupling equations in 
Lorentzian approximation read (note that z = 3/2) 



FIG. 4. Scaling function for the real part, Re[£'{(jj)], and 
imaginary part divided by the scaled frequency, Im [i'{il>)] /il>, 
of the generalized surface tension. 

Let us first discuss the vertex corrections resulting from 
the imaginary part of the generalized surfac e tens ion. Us- 

one can 



ing the mode coupling results from Section [II C , 
calculate the real and imaginary part of j/(k, w), as shown 
in Fig. |. 



^lor 



dy- 



1 



3/2 

y+ - 



y 



3/2 



(3.69) 



where y± — ^±y. This gives Vy^^ ~ 1.955. Next we take 
into account vertex corrections from the two-loop dia- 
grams. We have seen in section [I C that the two-loop 
contributions to v{q, 0) can be split into a propagator 
renormalization and a vertex correction. The former is 
already taken into account in the self-consistency scheme 
of mode coupling theory. Hence, in Lorentzian approxi- 
mation one may extend the mode-coupling approach in 
the following way 
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FIG. 5. Real part of the vertex correction versus jl, 
Im \y{jjL — ^) — {'{fi + ^)] /lo. The value of the scaling vari- 
able u) is indicated in the graph. 

Therefrom one deduces the real part of the vertex cor- 
rections 

Re [VkV^(K+; K) |k=o] + Aq = (3.68) 
-Aq/77i[j)(/i-|)-i)(/i + f)] /u , 



2 



l + V^^\p/q) 
p v{<l+)<l\ + v{q-)q^_ 



(3.70) 



with 



2y-l 



^lor 



1 

2 


+ y + X 


- 3 3 ■ 

yl+yl 




' — — 3 " 

yl+ymx\^ 



(3.71) 



where y± = |i ± ?/ ± a;| and y± = |^ ± ?/|. As can be 
inferred from Fig. ^ the vertex corrections V^^^ are not 
at all small as compared to the bare term 1. They vary 
from about —100% to -f-100% as a function of the ratio 
of the external momenta q and p. 

However, by inserting this vertex correction in the ex- 
tended mode coupling equation, Eq. ( 3.7C ), one finds 
^\o-[ ~ 2.044, which is merely less than 5% larger than 
the value obtained with the bare vertex. In addition, 
the correction tends to increase the amplitude ratio cal- 
culated in section HI C towards the value obtained from 
numerical simulations. This result clearly explains why 
mode coupling theory has been so successful in calculat- 
ing scaling functions for the noisy Bugers equation. 
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FIG. 6. Vertex correction V^^\y) from two-loop diagrams 
as a function of ln|j/|: dashed curve: y < 0, dot-dashed 
curve: y > 0. The sohd hne represents the average 

v^'Hy)=[v^'\y) + v^'\-y)] 12. 

It still remains a puzzle to us, however, that the vertex 
correction itself constitutes such a large correction to the 
bare vertex of the order of ±100%, while its effect on the 
correlation function is much less. Future investigations 
may concentrate on the extension of the two-loop vertex 
correction to higher orders, e.g. by including all ladder 
diagrams to the three-point vertex function, which may 
give us further insight in the validity and accuracy of the 
fairly simple mode-coupling approach. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have investigated the noisy Burg- 
ers equation using dynamical renormalization group and 
mode-coupling techniques. The renormalization group 
results show that there appears no singular two-loop con- 
tribution to the height-height correlation function. Upon 
including the two-loop vertex corrections into the mode- 
coupling approach we were able to show that their effect 
on the result for the correlation function is approximately 
5%, whereas the vertex corrections themselves are quite 
large. We suppose that this overestimates the actual 
vertex correction. In order to go beyond the two-loop 
vertex corrections one should possibly use an additional 
suitable resummation of the vertex correction, e.g., write 
down a "Bethe-Salpeter" type of equation for the ver- 
tices representing all ladder diagrams contributing to the 
three-point vertex functions. We leave this problem, and 
conceivable further extensions, for future investigations. 

As a complementary estimate for the vertex corrections 
we have used Ward identities to derive a relation between 
certain derivatives of the three-point vertex function and 
the renormalized noise amplitude and surface tension. 



Again, one finds that the effect of the vertex corrections 
may be of the order of a few percent. 

Based on the above estimates for the vertex correc- 
tions, we suppose that mode coupling theory yields very 
accurate results for the scaling functions, at least for the 
(l-t-l)-dimensional KPZ (noisy Burgers) equation. This 
conclusion is supported by the close agreement of the 
mode coupling results with those from numerical simula- 
tions of finite size systems. It remains unclear, however, 
whether the mode-coupling approach works as well in 
the general (d-l- l)-dimensional case, where the two-loop 
perturbation theory corrections do contain singular con- 
tributions. 
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APPENDIX A: GENERALIZED FLUCTUATION 
DISSIPATION RELATION 

In this appendix we collect several fluctuation dissi- 
pation theorems (FDT's), which are of importance for 
the dynamics of systems described by nonlinear Langevin 
equations. 

Let T be the time reversal operation: t —t. Then 
detailed balance (time inversion symmetry) implies [4^ ] 
that 



r exp [J-/^ ^Ft,]= exp [JZ*- ~ F^t,] 



(Al) 



where Ft — F[h{t)] is the stationary probability distri- 
bution function and 



Jl^Kh] = / dx 



dt 



h{x,t)Dh{x,t) 



with 



(A2) 



(A3) 



Without loss of generality, one can assume that the field 
h(x,t) is even or odd under time reversal, that is 



Th = eh, £ = ±1 



(A4) 



Here h is odd under time reversal. The stationary dis- 
tribution Pst[^] = e~^[''l is characterized by the "free 
energy" 
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The time reversal symmetry implies that 

Sh{-t) 



Th(t) = -e h{-t) 



(A5) 



(A6) 



Now one uses the causality property of the response 
functions, {h{ti)...h{tk)h{ii)...h(tk)) = 0, if one ij) allti. 
Then for example 



{h{t)h{0)) = for t < 



(A7) 



With the time reversal operation and Eq. (A6) it follows 
then from Eq. (|A7|) for t < that 







Upon redefining t = —t, one obtains for t > 

(Mt)M0)) = ew(/.wM). 



(A8) 



(A9) 



The same arguments can be repeated for 
{h{ti)...h{tk)h{ii)) with ii > aU tj. The result is 



{h{h)...h{tk)h{h)) 



eih,{t,}){hih)...h{tk) 



SFlhjh)] 
5h{ti) 



(AlO) 



where <d{ti,{tj}) is an obvious generalization of the 
function. Note that these generalized FDT's are for the 
cumulants and not for the vertex functions. In particular 
we get 



Gii(fc,t) = e(t) — Go2(fc,t). 



(All) 



Further identities can be written down in a completely 
analogous way Pq,E3l . 



APPENDIX B: TWO-LOOP PERTURBATION 
THEORY FOR THE TWO-POINT VERTEX 
FUNCTIONS 

This appendix comprises the Feynman diagrams to 
two-loop order for the (l+l)-dimensional Kardar- 
Parisi-Zhang equation, and the corresponding momen- 
tum integrals. The integrations over the internal fre- 
quencies have already been performed using the residue 
theorem. 



We start with a list of the contributions to two-loop 
order to the fully wavevector- and frequency-dependent 
two-point vertex function r^^(g,w). The other non- 
vanishing vertex function F,- r (q, oj) can be calculated in 
a similar fashion (see Ref. B), or sim ply be obtained 
via the fluctuation-dissipation theorem ( p.8[ ). In writing 
down the diagrammatic expansion for the dynamic func- 
tional one has to take into account restrictions which 
follow from causality. In section II we have not explicitly 
included the Jacobian J[h] = 'D[r]]/'D[h], which depends 
on the discretization of the Langevin equation (needed 
to give a proper definition to the path integral). As can 
be shown quite generally |^^, the Jacobian cancels the 
equal-time contractions of the field h and the response 
field h. Keeping this in mind (or by choosing a discretiza- 
tion with the Jacobian equal to 1), one can omit the Jaco- 
bian in the dynamic functional. The Feynman diagrams, 
which account for the restrictions imposed by causality, 
are depicted in Fig. |7[ 




FIG. 7. Feynman diagrams for the two-point vertex func- 
tion r,-^(q,tj) to two-loop order. 



Introducing the abbreviations q± 
expressions read: 



{q/2) ± p, q± — (g-/2) ± k, and q± — q± k, the corresponding analytical 



(a) -t- (6) — iuj + v^q^ + 



2vo ' 



p iuj + i^oq+ + i^oQ- 



(2.1) 
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(c) = 
(e) = 

(/) 
{9)- 

(h) 

ii) 



1 



Q+QZ 



jdk 



1 



1 1 



1 



q+q- 



q+qt 



ik 



2*^0 Jp + ^oq+ + ^oq- Jk ^oq- + '^oq+ + ^oq- V ^o9- ' «w + uoq\ + uoq\ + i/Q?^ 



+ - 



9+ 



2mi- 



X 1 + 



q- 



q- f q+q- 



q+k 



"0 Jp 
^0 Jp 



iu + i/oq+ + voqt Jk [«w + ^oq+ + I'oa-] [«w + voq\ + Voq^_ + i/ofc^] 



9+ 



iu> + voq% + voq^_ J^. [uoq'i + V(y$_ + fofc^] [iuj + i^og^ + uoqi] 



X ( 1 + - 

UjJ ■ 



+ yoq+ + t^oq- + yok"^ J 
^"Mq f q q_ f l± 



(2.2) 
(2.3) 

(2.4) 
(2.5) 
(2.6) 

(2.7) 
(2.8) 



(2.9) 
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